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Achieving coherent population transfer in the solid-state is challenging compared to atomic systems due to 
closely spaced electronic states and fast decoherence. Here, within an atomistic pseudopotential theory, we 
offer recipes for the stimulated Raman adiabatic passage in embedded silicon and germanium nanocrystals. 
The transfer efficiency spectra displays characteristic Fano resonances. By exploiting the Stark effect, we 
predict that transfer can be switched off with a dc voltage. As the population transfer is highly sensitive to 
structural variations, with a choice of a sufficiently small two-photon detuning bandwidth, it can be harnessed 
for addressing individual nanocrystals within an ensemble. 
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The control of the dynamics of quantum systems us- 
ing coherent optical beams lies at the heart of quantum 
information technologies Among several alternatives, 
the stimulated Raman adiabatic passage (STIRAP) of- 
fers a certain degree of robustness in atomic systems with 
respect to laser parameter fluctuations & Its solid-state 
implementation was recently achieved in Pr +3 :Y2SiC>5 
crystal,—^ and Tm +3 :YAG crystal,— all at cryogenic tem- 
peratures. Next obvious milestone is to demonstrate STI- 
RAP in nanocrystals (NCs) embedded in a host lattice. 
Compared to rare-earth doped ions in inorganic solids,— — 
NCs bring further flexibility in the design of the func- 
tional units by tailoring the physical parameters such as 
material composition, size, shape, and strain, together 
with the external fields. However, the challenge with NCs 
is that the charge degrees of freedom is more susceptible 
to decoherence compared to the atomic or "impurity" 
systems 

In this work we consider silicon and germanium NCs 
embedded in silica. Our aim here is to explore from a 
theoretical standpoint the feasibility as well as the intri- 
cacies of STIRAP in this system. To set the stage, first 
we need to address the constraints imposed by decoher- 
ence on our system. The ultimate decoherence mecha- 
nism is the radiative recombination. The typical radia- 
tive lifetimes of Si and Ge NCs are in the microsecond 
range whereas for direct band gap semiconductors this is 
in the nanosecond ranged Another recombination chan- 
nel, in case multiple electrons get excited by a strong 
laser pulse, is the Auger process. According to our re- 
cent theoretical estimation for the excited-electron con- 
figuration of Auger recombination in Si and Ge NCs, this 
lifetime is in the range of sub-nanoseconds.— An even 
more critical decoherence channel in NCs is the acoustic 
phonon scattering!^ For the case of InGaAs quantum 
dots, Borri et al. have demonstrated close to radiative 
limit linewidth at 7 K, corresponding to a dephasing time 
of 630 ps»ii In Si NCs, Sychungov et al. have shown that 
the linewidth can also be as sharp as direct band gap 
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materials, reaching 2 meV at 35 A similar system is 
the excited Rydbcrg states^ of phosphorus-doped silicon 
having a spatial extend of ~10 nm for which a dephas- 
ing time of ~320 ps is very recently predicted.— Guided 
by these reports, we aim for a complete STIRAP in less 
than 300 ps so that at sufficiently low temperatures of 
a few Kelvins this can beat the decoherence clock in Si 
and Ge NCs. Admittedly, this is a cautiously optimistic 
estimate, nevertheless a worse case can still be accom- 
modated by further scaling the pulse widths and laser 
powers accordingly; thanks to high-field tolerance of sil- 
ica embedded NCs. 

Our theoretical model involves the atomistic descrip- 
tion of the system within a supercell geometry of sev- 
eral thousand atoms most of which are the surround- 
ing matrix atoms. Initially, nearly spherical NCs in 
C% v point symmetry are considered, and in the final 
part consequences of shape deformation are discussed. 
The local potential is represented as a superposition 
of screened semiempirical pseudopotentials of the con- 
stituent atoms;— the spin-orbit interaction is particularly 
included, since this coupling among closely spaced levels 
can potentially affect the selection rules and hence the 
transfer efficiency. A dc electric field is also accounted 
nonperturbatively for the Stark field analysis. The ex- 
citonic effects are ignored as the confinement energy 
dominates for small NCs4 To solve the single-particle 
Schrodinger equation for such a large number of atoms 
with sufficient accuracy up to the highly excited states, 
we make use of the linear combination of bulk bands 
approach^ 

The population transfer is built on this atomistic 
electronic structure as shown in the insets of Fig. 1. 
The electric dipole coupling is used for the interaction 
with the pump and Stokes beams. Unlike atomic sys- 
tems, in the case of NCs we have to consider multiple 
intermediat o 16 i 17 and final states. On the other hand we 
assume a single initial state, namely, the highest occu- 
pied molecular orbital (HOMO). However, by selecting 
the interaction parameters accordingly (see Table I) we 
have assured that the maximum probability of finding 
a sub-HOMO electron in the conduction band is quite 
negligible. Our computations show that had the interac- 
tion individually involved any such sub-HOMO electron 
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as the initial state, its transfer probability to the conduc- 
tion band would have an upper bound of 10~ 9 %. Finally, 
since the intermediate state is not populated for the case 
of an ideal STIRAP, we neglect many-body effects due 
to Pauli blocking. 

We base our discussions on a spherical 2.1 nm diameter 
Si NC and a 1.5 nm Ge NC, purposely selected because 
of their close band gaps of 2.74 eV and 2.80 eV, respec- 
tively. Considering the oscillator strength of the transi- 
tions, we adopt different schemes for the two NCs (see 
right insets in Fig.[T]). For 1.5 nm Ge NC the HOMO to 
lowest unoccupied molecular orbital (LUMO) transition 
is quite strong, therefore, we utilize the ladder scheme 
where the electron is transferred from the HOMO state 
to the LUMO+28 state (0.837 eV above LUMO) via the 
LUMO state. For 2.1 nm Si NC, we utilize a A scheme 
where the electron is transferred from the HOMO state 
to the LUMO state via the LUMO+30 state (0.937 eV 
above LUMO). For the pump and Stokes pulses, we use 
counterintuitivcly time-ordered Gaussian profiles^ giving 
rise to equal peak Rabi frequencies with values 0.35 and 
0.55 THz for Si and Ge NCs, respectively. Both pulses 
are linearly polarized but in certain directions chosen to 
optimize the transfer. We refer to Table I for the other 
laser parameters. 

TABLE I. Laser parameters optimized for STIRAP for the 
2.1 nm Si, and 1.5 nm Ge NCs. The incident electric fields are 
specified for free-space medium. Delay refers to time between 
the peaks of the Stokes and pump pulses both with Gausian 
profiles. 
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The transfer efficiency is plotted in Fig. [TJ In all the 
figures, detuning refers to that of the pump pulse; in 
the case of two-photon resonance, the Stokes pulse is as- 
sumed to be detuned by the same amount from its tar- 
geted transition. In the other case, the Stokes pulse is 
kept resonant. The first observation we make is that 
the transfer efficiency is superior when the two-photon 
resonance condition is satisfied. This fact is well docu- 
mented in literature.— When the system has two-photon 
detuning, the central peak shrinks down into a plateau, 
where the transfer efficiency is essentially constant, and 
suddenly drops to zero (see upper inset). 

Another noteworthy aspect of Fig. Q] is the presence of 
side peaks. The detuning values for these peaks coincide 
with the resonance condition restored with a neighbor- 
ing intermediate state, marked by the blue vertical bars 
in Fig. 1. As the laser parameters were optimized for 
the central peak and not for the neighbors, these peaks 
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FIG. 1. (Color online) The population transfer efficiency as 
the pump laser is detuned for the 1.5 nm Ge NC (top) and 
the 2.1 nm Si NC (bottom), with (8 = 0) and without two- 
photon resonance (<5 7^ 0). The left inset in the upper graph 
is a close-up for the central peak. The vertical blue lines 
on the detuning axis show the energies of the intermediate 
states. Insets on the right show the electronic states and laser 
energies. 



are usually not as tall or wide. Also since there are no 
states immediately below LUMO (ladder scheme) , for the 
1.5 nm Ge NC there are no such peaks for negative de- 
tuning. The peaks display the asymmetrical well-known 
Fano lineshapei 18 ' 19 In a similar context, this was ob- 
served in the tunelling-induccd transparency in quantum 
well intersubband transitions^ It arises from two paths 
interfering with opposite phase on one of the two sides of 
the resonance. This is illustrated in the lower left inset of 
Fig. [1] For small detunings the calculations made with- 
out considering any neighbors agree very well with the 
full calculation. However, as soon as there is enough de- 
tuning to transfer the electron through one of the neigh- 
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FIG. 2. (Color online) Effect of dc Stark field on the overall 
transfer efficiency and maximum intermediate state popula- 
tion for the 2.1 nm Si NC. Lines are only to guide the eyes. 



boring states, the neighbors-removed treatment cannot 
reproduce the dip right before the second peak in the 
solid line which occurs due to the interference between 
the chosen intermediate state and its neighbor. 

Next, in Fig. [2J we investigate the effect of an exter- 
nal dc electric field (see the inset). Note that we quote 
the matrix dc electric fields, viz., outside the NC in the 
embedding region of silica. Here, for the 2.1 nm Si NC, 
we observe that for small fields, the transfer efficiency 
is not affected, but the intermediate population pile-up 
increases. After a critical field of 0.35 MV/cm, transfer 
efficiency rapidly drops to zero. To identify its origin, we 
first assured that the system is robust against changes in 
the dipole matrix elements, hence Rabi frequencies are 
not significantly altered by the dc field. On the other 
hand, NC energy levels undergo significant Stark shifts, 
the valence states being more so compared to conduction 
states, as revealed by our recent work<2I We checked that 
in this field range it does not give rise to a level cross- 
ing between HOMO and the lower-lying states. Hence, 
the primary mechanism responsible for this switching is 
the Stark shift-induced two-photon detuning, something 
STIRAP is very sensitive to. For the zero-field case, a 
detuning of 0.3 meV is enough to destroy STIRAP. This 
value of Stark shift is reached at 0.5 MV/cm, after which 
the population transfer is quenched. 

The highly critical two-photon bandwidth is control- 
lable by the time delay between Stokes and pump pulses 
which is illustrated in Fig. |3l As the overlap between 
these pulses is reduced the two-photon bandwidth first 
increases up to a delay of 80 ps, beyond which it re- 
tracts back, as expected from the fundamental principles 
of STIRAP that demands a non-zero overlap between the 
two pulses<£ The upper inset shows the build-up of the 
intermediate-state population away from the two-photon 
resonance. 
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FIG. 3. (Color online) Two-photon detuning versus transfer 
efficiency for different time delays for the 2.1 nm Si NC. Upper 
inset: the maximum probability of finding the electron in any 
of the intermediate states throughout the transfer process. 
Lower inset: the effect of NC ellipticity on the two- photon 
detuning and transfer efficiency (Lines are guides to the eyes); 
the horizontal dashed line marks the critical 0.3 meV two- 
photon detuning. 

Finally, we focus on the NC's structural sensitivity. 
Starting with its shape, we consider NCs of the same 
number of atoms (hence, same volume under zero strain) 
but with different asphericities as quantified by the el- 
lipticity parameter, e. Denoting a and b as the equato- 
rial radii and c as the polar radius, for oblate spheroids 
(a = b > c) e = \J\ — c 2 /a 2 , whereas for prolate 
spheroids (a = b < c) we define it to be negative as 
e = —y/l — a 2 /c 2 . The lower inset of Fig. [3] vividly dis- 
plays the fact that the transfer is lost as the shape of the 
NC is deformed from the originally targeted geometry 
to which STIRAP was optimized (here spherical). The 
repositioning of only two NC surface atoms (e=0.3 case) 
is enough to displace the electronic states away from the 
tolerable two-photon detuning window. Likewise, we ob- 
served that an incremental change in the size of the NC 
by including the next shell of atoms (not shown) results 
in a similar loss of transfer. These indicate that prac- 
tically the intended STIRAP will be locked only to the 
single NC that it is tuned to. 

In conclusion, we provide a theoretical insight for STI- 
RAP in small Si and Ge NCs. Due to dense electronic 
states it displays a train of Fano resonances. The trans- 
fer can be abruptly switched off with a dc voltage by 
introducing Stark shift that sufficiently detunes the two- 
photon resonance. Finally, we demonstrate the sensitiv- 
ity of the transfer efficiency with respect to the structure 
of the NC which can be instrumental in addressing a sin- 
gle NC among an ensemble having inherent size, shape 
and even local strain fluctuations. 
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